IMF isotopic properties in semi-peripheral collisions at Fermi energies 
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We study the neutron and proton dynamical behavior along the fragmentation path in semi- 
peripheral collisions: 58 Fe+ 58 Fe (charge asymmetric, N/Z = 1.23) and 58 Ni+ 58 Ni (charge symmetric, 
N/Z = 1.07), at 47 AMeV. We observe that isospin dynamics processes take place also in the 
charge-symmetric system 58 Ni+ 58 Ni, that may produce more asymmetric fragments. A neutron 
enrichment of the neck fragments is observed, resulting from the interplay between pre-equilibrium 
emission and the phenomenon of isospin-migration. Both effects depend on the EoS (Equation of 
State) symmetry term. This point is illustrated by comparing the results obtained with two different 
choices of the symmetry energy density dependence. New correlation observables are suggested, to 
study the reaction mechanism and the isospin dynamics. 

PACS numbers: 25.70-7, 21.30.Fe, 24.10.Cn, 21.65. +f 
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Collisions between heavy ions with large isospin asym- 
metries, made possible by the recent radioactive beam 
developments, represent a very efficient way to probe the 
structure of nuclear- EoS (Equation of State) symmetry 
term. In particular the symmetry energy behavior is in- 
fluencing very dissipative reaction mechanisms, such as 
fragmentation processes, leading to important effects on 
fragment composition. 

In central heavy ion collisions at intermediate ener- 
gies the spinodal decomposition has been proposed as a 
possible mechanism for fragment formation, see Ref.Q 
for a recent review. According to this description, frag- 
ments should reflect the properties of the low density 
phase, where they are formed. In charge asymmetric sys- 
tems, the isospin distillation, i.e. the formation of more 
symmetric fragments surrounded by a neutron richer di- 
lute phase, takes place. Here we will focus on fragmen- 
tation in semi-peripheral collisions, where intermediate 
mass fragments (IMF) are mostly produced in the over- 
lap zone (the neck region) between projectile-like and 
target-like fragments (PLF /TLF, the spectator region) , 
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We will discuss the influence of the symmetry energy on 
the features of this fragmentation mechanism. The pres- 
ence of a density gradient between the neck (low density) 
and the spectator (high density) regions affects the N/Z 
of fragments in a different way with respect to the spin- 
odal decomposition mechanism [2(|, l2ll |22| . Moreover 
since the isoscalar density gradients are ruling the isospin 
transfer through the density dependence of the symmetry 
energy, we see that measurements of isospin observables 
in semiperipheral collisions will directly probe the slope 
of the symmetry term around saturation, of large impor- 
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tance for the structure of neutron-rich nuclei [23| . 

We consider the reactions 58 Fe + 58 Fe (charge asym- 
metric N/Z = 1.23) and 58 Ni + 58 Ni (charge symmetric 
N/Z — 1.07). In fact in these collisions, due to the uni- 
form N/Z distributions, we do not have isospin gradients 
initially. We show that the study of the full reaction dy- 
namics and of the possible occurrence of density gradients 
is essential in order to understand the isospin dynamics. 
We will finally suggest the measurement of some correla- 
tions between IMF properties, like mass, isospin content 
and "alignement" , particularly sensitive to the reaction 
mechanism and the isospin transport. 

To get an insight into the behavior of neutrons and pro- 
tons in asymmetric matter one can consider the density 
dependence of neutron and proton chemical potentials: 
p q = de(p q , p q >)/dpq, q = n,p,e being the energy density. 
We recall that this quantity contains all contributions to 
the energy per particle (kinetic, potential and simmetry 
energy). In Fig ^ we report the density dependence of 
the n,p chemical potentials below normal density, where 
we expect that fragment formation takes place 7], for a 
system with asymmetry / = (N — Z)/A = 0.2 and for 
two choices of the iso- EoS. We refer to an asy — stiff 
EoS when we consider a potential symmetry term that 
is linearly increasing with nuclear density. We refer to 
an asy — soft EoS when the potential symmetry term 
increases up to a saturation around normal density, and 
then eventually decreases [2(3, El • 

Since particles move towards the minimum of the 
chemical potential, it is possible to observe that there 
exists a density window, roughly po/2 < p < po from 
Fig^ where neutrons and protons can move in opposite 
directions: protons move towards a higher density region 
while neutrons move towards a lower density region. This 
phenomenon, called isospin migration |2l| , causes a neck 
neutron enrichment, due to the density gradient between 
the low density neck region and the spectator matter. We 



2 



y 1 1 1 I 1 1 1 1 I 1 1 

■ \ \ 
- u \ 

A \ s " 


,,,„,,_ 


\ \ 
V \ 

\ \ 
\ \ 
\ \ 
\ \ 






,-,-f , , , , 1 , ■ 



0.10 0.15 
-3s 



p(fm- 3 ) 

FIG. 1: Density dependence of the chemical potential for neu- 
trons (upper curves) and protons (lower curves) for an asy- 
stiff (solid lines) and asy-soft (dashed lines) EoS with asym- 
metry parameter / = 0.2. 
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FIG. 3: Asymmetry vs. charge of each fragment arising from 
the simulation of the reaction 58 Fe + 58 Fe (a) and 58 Ni + 
58 Ni (b) with an asy — stiff EoS. Horizontal dashed lines 
are the initial asymmetries of the colliding systems. 
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FIG. 2: Density contour plots on the reaction plane for two 
different events of the reaction 58 Fe + 58 Fe at 47 AMeV, 
b r = 0.5. The box size is 40/m. The numbers inside the box 
represent the time-steps in fm/c. 



stress that this mechanism is different from the recently 
investigated isospin diffusion and N/Z equilibration pro- 
cesses in peripheral collisions [3 E3- that instead are 
due essentially to the presence of isospin gradients. 

Reactions have been simulated by considering a 
stochastic extension of the microscopic transport equa- 
tion BNV (Boltzmann-Nordheim- Vlasov), following a 
test-particle evolution on a lattice [2(| I27L |2S| . We con- 
sider a beam energy of 47 AMeV, and reduced impact 
parameter b r = b/bmax — 0.5 (semiperipheral) . Using 
the asy — stiff EoS 40% of the events produce at least 
one IMF in the neck region (ternary events). In Fig. |2 
density contour plots on the reaction plane are presented 
for two different events coming from the reaction 58 Fe + 
58 Fe (b r = 0.5, 47 AMeV). 



The first row shows an event in which a fragment forms 
with a time delay and in a space region correlated to a 
target-like nucleus, unlike the event in the second row, 
which shows a rapid fragment formation in a region that 
is not correlated to any spectator remnants. Thus frag- 
ments are formed according to a variety of mechanisms 
and we will see how this influences their properties. 

Fig. |2| reports the N/Z ratio of each fragment vs. the 
charge Z at the freeze-out time, obtained at the interme- 
diate impact parameter b r = 0.5. 

The residual PLF/TLF nuclei (large Z range) show a 
different behavior in the two reactions: we note, in fact, 
that the points of Fe system are along the dashed line, 
that represents the initial system asymmetry, while for 
Ni reaction points lie above that line. The IMF (low 
Z range) behavior is similar for the two reactions: the 
points, for both reactions, lie above the dashed line, al- 
though for the reaction 58 Ni + 58 Ni the difference be- 
tween the N/Z of IMF's and large PLT/TLF residues 
seems to be less pronounced. 

For a better understanding of Fig|3 we have to con- 
sider what happens during the pre-equilibrium phase 
[1 12^. l3fH l3l| . The asymmetry of the di- nuclear neutron- 
rich system changes from 1.23 (initial value) to 1.22 (at 
t = 100 fm/c, instant in which fragments start to form) 
since 14 neutrons and 11 protons are evaporated, while 
the di- nuclear neutron-poor system changes from 1.07 to 
1.12 as a consequence of a larger proton evaporation, due 
to the Coulomb repulsion (it loses 13 protons against 12 
neutrons), becoming an asymmetric system. We can con- 
clude that: 

i) In the neutron rich reaction, the neutron emis- 
sion due to pre-equilibrium goes in the same direction 
of the neck neutron enrichment, caused by the isospin- 
migration. So finally we observe slightly more symmet- 
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TABLE I: Asymmetry evolution of the residual nuclei arising 
from binary and ternary events. 
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t = 


t = 100/m/c 


t = 200/m/c 


58p e + 58 Fe 


1.23 


1.22 


1.23 binary 








1.19 ternary 


5S Ni + 58 Ni 


1.07 


1.12 


1.17 binary 








1.125 ternary 
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FIG. 4: Average asymmetry vs. Z for nuclei from the reac- 
tions 58 Fe + 58 Fe (left panel) and 58 Ni + 58 Ni (right panel) for 
an asy — stiff (circles) and asy — soft (triangles) EoS. Hor- 
izontal dashed lines represent initial asymmetry of colliding 
systems [2l| . 



FIG. 5: Contour plots of the parallel velocity-mass distribu- 
tion for the two systems: 58 Fe + 58 Fe (left) and 58 Ni + 58 Ni 
(right). Arrows indicate the positions of average PLF and 
TLF parallel velocity. 
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ric PLF /TLF residues accompanied by neutron-richer 
IMF's; 

ii) In the neutron-poor collision, the larger pre- 
equilibrium proton emission enhances the N/Z of the 
dinuclear composite system. The acquired asymmetry 
is then transferred to the neck region. 

To better disantangle between the two mechanisms, 
one can study the charge composition of residues distin- 
guishing between binary and ternary events (Table [IJ. 

For both reactions we note that the N/Z ratio of resid- 
ual PLF '/TLF nuclei in ternary events is lower than 
the value for binary events, since in the latter case the 
isospin-migration effect does not apply. The isospin dy- 
namics effect is rather evident from the comparison with 
the asymmetry values at the end of the pre-equilibrium 
phase (t = 100/m/c in the Table). For the Fe + Fe 
system the N/Z of residues changes from 1.22 to 1.19, 
in ternary events; for the Ni + Ni reaction this differ- 
ence is not so evident (from 1.12 to 1.125) because the 
isospin-migration competes with proton evaporation. On 
the other hand, in binary events, we note the neutron en- 
richment of residues in the Ni reaction, due to a favorite 
proton pre-equilibrium emission. 

In Fig. 0| the average asymmetry of products arising 
from the two reactions is shown for the two choices of 
the symmetry energy parameterization, asy — stif f and 
asy — soft EoS. 

We note that fragments are more symmetric in the asy- 
soft case. This difference can be explained with the dif- 



FIG. 6: Contour plots of the N/Z distribution versus the emis- 
sion angle. Asy — stiff choice. 



ferent behavior of the chemical potential in the two EoS, 
see Fig. ^ that will induce a different isospin-migration 
effect. For an asy-soft EoS, the proton chemical poten- 
tial varies not so much in the region where the fragments 
form (from 0.08 to 0.15 fm~ 3 ), while for neutrons there 
is a significant slope (however smaller than in the asy- 
stiff case) . Therefore neutron enrichment of the neck will 
still cause an increase of the fragment asymmetry, but 
this increase will be smaller than in the asy-stiff case, 
where protons can even migrate out of the neck region. 
The fact that the N/Z of large fragments is smaller in 
the asy-soft case is due to pre-equilibrium effects, since 
with the asy-soft EoS more neutrons are emitted due 
to the more repulsive mean field below normal density, 
[S3, HIE! EI- On the other hand, the larger difference 
between IMF and spectator fragment asymmetries with 
the asy-stiff parametrization is a consequence of a more 
effective isospin migration, see the recent |33| . 

It is interesting to look also for some correlations be- 
tween asymmetry, mass, velocity and direction of the out- 
going fragments. As we will show, in this way we can even 
disentangle between features just related to the reaction 
mechanism and observables more sensitive to the isospin 
dynamics and symmetry energy. 
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FIG. 7: Average IMF mass as a function of the distance from 
the PLF-TLF axis at the freeze-out time for 47 AMeV, b r = 
0.5 collisions. Left panel: Fe + Fe. Right panel: Ni + Ni. 
Empty squares: asy — soft symmetry term. Full squares: 
asy - stiff. 
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FIG. 8: Ratio of the IMF yields, with N/Z larger and smaller 
than the value a obtained just after pre-equilibrium emission, 
as a function of the distance from the PLF — TLF axis at the 
freeze-out time for 47 AMeV, b r = 0.5 collisions. Left panel: 
Fe + Fe. Right panel: Ni + Ni. Empty squares: asy — stiff 
symmetry term. Full squares: asy — soft. 



In Fig|S] we present the behavior of the velocity along 
the beam direction versus the IMF mass. Lighter frag- 
ments are emitted at all angles and they have larger ve- 
locities with respect to more massive IMF's that are 
more correlated to the spectator matter and are emitted 
on longer time scales. 

The analysis of the N/Z distribution versus the emis- 
sion angle (see FigEJl reveals that larger fluctuations are 
present close to forward and backward angles. This indi- 
cates that IMF's more correlated to the spectator mat- 
ter may become more neutron-rich, since they interact for 
a longer time with the system and the isospin transport 
mechanism becomes more effective. 

The presence of correlations between mass, neutron 
content and kinematical observables can be better evi- 
denced by studying IMF properties vs. the "alignement" , 
i.e. as a function of the distance d from the PLF — TLF 
axis at the freeze-out time. 

In Fig[7|we plot the behavior of the IMF average mass 
vs. the distance d for the n-rich (FeFe, left panel) and 



the n-poor (NiNi, right panel) collision. The calcula- 
tion is performed with the two choices of the density de- 
pendence of the symmetry term. In both cases we see 
a clear increase of the IMF masses with the fragment 
"alignement", independent of the stiffness of the used 
I so — EoS. This appears a rather general feature of the 
neck — fragmentation mechanism: light fragments are 
emitted at earlier times and are not much driven by the 
spectator residues (PLF /TLF) Q. This is also in line 
with the fact the light fragments may reach larger veloc- 
ities (see FigjSJ. 

As a consequence we would expect a different amount 
of isospin migration vs. the IMF alignement, and now 
such correlation should be I so — Eos dependent. In Fig[5] 
we report the ratio of the IMF yields with the N/Z 
larger and smaller than the value, a, reached just after 
pre-equilibrium emission (see Table I), plotted vs. the 
alignement distance d at freeze out. For both systems we 
have a nice increase of the neutron enrichment with the 
alignement. As expected from the previous discussion on 
the physics of the isospin migration, the effect is more 
evident in the asy — stiff choice. 

We stress again that neck — IMF's always present a 
neutron enrichment, even in the case of a n-poor sys- 
tem. The latter paradox is due to the pre-equilibrium 
isospin dynamics. In fact, as anticipated above, due to 
pre-equilibrium emission, the system will loose some pro- 
tons and acquire a N/Z larger than the initial one. Then, 
before the di-nuclear system reseparates, the neutron ex- 
cess is transferred to the neck region that is at lower 
density. Both effects, fast proton emission and neutron 
transfer, are connected to the symmetry term of the EoS. 
Some evidence has been found in recent data on 58 Ni in- 
duced fragmentation 0, 0, . 

In conclusion, we have shown that isospin dynamic pro- 
cesses appear even in systems with initially uniform spa- 
tial asymmetry distribution, such as 58 Fe + 58 Fe and 58 Ni 
_l_ 58 Ni. From a chemical point of view, we do not expect 
N/Z gradients which can induce asymmetry variations. 
These variations are instead caused by density gradients 
during the reaction dynamics since the symmetry term 
of the EoS introduces a different behavior of the chem- 
ical potentials for neutrons and protons with respect to 
density. So, when the collision happens, the spatial dis- 
tribution of the isoscalar density will induce variations 
even in the isovector density. 

We have revealed an interesting correlation, typical of 
the neck — fragmentation mechanism, between IMF 
masses and corresponding alignement to the axis joining 
the PLF — TLF residues. This suggests a time-hierarchy 
in the mid-rapidity fragment production with the lighter 
clusters formed at earlier times. When we combine to the 
isospin migration dynamics, an I so — Eos sensitive ob- 
servable results to be the correlation between the neutron 
excess of IMF's and the relative alignement. 

The reactions 58 Fe + 58 Fe and 58 Ni + 58 Ni have been 
studied experimentally, focussing on central collisions 
and adopting a statistical description of multifragmen- 
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tation in [34ll35|. or in semi-peripheral collisions in [l5j . 
Asymmetry effects, as the ones described here, have been 
noted in (T^J, even in the second system. In this paper 
we suggest new interesting correlation analyses that can 
be pursued selecting semi-peripheral events and that ap- 



pear particularly appropriate in order to investigate the 
isovector structure of the EoS as well as to shed some 
lights on the mechanism responsible for fragment pro- 
duction. 
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